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INFLUENCE OF RANDOM BULK INHOMOGENEITIES ON QUASI-OPTICAL 

CAVITY RESONATOR SPECTRUM 
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Institute for Radiophysics and Electronics, 
National Academy of Sciences of Ukraine, 
12 Proscura St., 61085 Kharkov, Ukraine 

We suggest the statistical spectral theory of oscillations in quasi-optical cavity resonator filled 
with random inhomogeneities. It is shown that inhomogeneities in the resonator result in inter- 
mode scattering leading to the shift and broadening of spectral lines. The shift and broadening 
of each line essentially depends on frequency distance to adjacent spectral lines. With increasing 
the distance the influence of inhomogeneities sharply reduces. The solitary spectral lines which 
have the distance to the nearest lines quite large is slightly changed due to small inhomogeneities. 
Owing to such selective influence of inhomogeneities on the spectral lines the effective spectrum 
rarefaction appears. Both the shift and broadening of spectral lines as well as spectrum rarefaction 
in quasi-optical cavity millimeter wave resonator were detected experimentally. We found out that 
inhomogeneities result in stochastization of the resonator spectrum in that mixed state appears, i.e. 
the spectrum acquires both regular and random parts. The active self-oscillator system based on 
the inhomogeneous quasi-optical cavity millimeter wave resonator with Gunn diode was studied as 
well. The inhomogeneous quasi-optical cavity millimeter wave resonator (passive and active) can 
serve as a model of semiconductor quantum billiard. Based on our results we suggest using such 
billiards with spectrum rarefied by random inhomogeneities as an active system of semiconductor 
laser. 

PACS numbers: 05.45.Mt, 42.25.Dd, 42.60.Da, 32.30.Bv 



INTRODUCTION 



The problem of electromagnetic wave propagation in 
random inhomogeneous media has been the issue of the 
day for several decades. Numerous publications are de- 
voted to the analysis of different aspects of this problem 
(see [1,2] and references therein). In this scientific area 
the subject of research is normally the scattering of elec- 
tromagnetic (acoustic) waves in unbounded or partially 
confined systems (for example, in waveguides) which con- 
tain random inhomogeneities. Researches on wave propa- 
gation in open statistically unregulated systems are stim- 
ulated by numerous applications for long-distance signal 
transmission in both radio and optical wave ranges. 

The problem of electromagnetic oscillations in con- 
fined systems with random inhomogeneities has signifi- 
cant meaning as well. At the same time, by now it has 
not been developed enough both from theoretical and 
experimental side. It is caused by, in particular, modern 
theoretical studies of electromagnetic wave propagation 
and scattering in random inhomogeneous media use as- 
sumption about the isotropy and homogeneity scatter- 
ing medium filled confined systems [1,2]. Therefore, they 
cannot be applied in general to the study of confined 
systems such as cavity resonators. The experimental 
study also requires new approach to study electromag- 
netic oscillations in inhomogeneous cavity resonators in 
wide wavelength band. 
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In our recent paper [3] we suggested a new spectral 
approach to study confined systems with random inho- 
mogeneous. There spectral properties of cavity spherical 
quasi-optical millimeter wave resonator filled with ran- 
dom sapphire particles were considered. The sapphire 
particles with dimensions of the order of operating wave- 
length affect significantly the resonator spectrum. The 
spherical frequency degeneration is completely removed 
and the spectral lines thus have chaotic distribution on 
the frequency scale. The lines became wider and the 
quality factor is correspondingly decreased. The analo- 
gous broadening of spectral lines caused by random in- 
homogeneities has been detected in Refs. [4,5], where 
resonators with random rough boundaries were studied. 

Recently the question about the influence of random 
inhomogeneities on the resonator spectrum has attracted 
a great attention in relation to the design of lasers 
on open micro-resonators [6]. In such resonators, the 
whispering-gallery modes with super-high quality factors 
can be excited. The possibility to realize such quality 
factors depends essentially on the number of inhomo- 
geneities (roughness level) on the resonator boundaries. 

It should be realized that cavity resonator with ex- 
tremely small dissipation loss is the almost Hermitian 
system whose oscillations are formed due to the restricted 
motion of electromagnetic waves. Therefore, the res- 
onator spectrum is discrete. Random inhomogeneities 
whether bulk or on the resonator boundaries not intro- 
ducing additional dissipative loss in general have not to 
broaden the spectral lines. Nevertheless, the broadening 
of lines in the spectrum effects and its stochastization 
have been detected in experiment [3]. 
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In this connection the question arises: what is the 
physical mechanism of spectral lines broadening and 
spectrum stochastization in the random inhomogeneous 
resonator? The study of this question is one of the rea- 
sons of the present paper. Another motivation is relevant 
to the nanoelectonics problem. Recently the great em- 
phasis is the study of a new type of nanoelectron systems 
(so-called zero-dimensional ones). The charge carrier mo- 
tion in these systems has space restriction in all three 
dimensions. The peculiar example of zero-dimensional 
system is a quantum dot (QD). QD is a semiconductor 
area by the size of order of 10 nm with electron (hole) 
conductivity and restricted from outer area by potential 
barrier. Owing to finite charge carrier motion the energy 
spectrum in QD is discrete and the number of spectral 
levels due to small QD size is relatively few. 

Recently an elegant method of QD array implementa- 
tion has been designed. This method is based on the self- 
organization effect in strained double GaAs heterostruc- 
tures [7-9]. The usage of QD array as an active medium 
permits to design lasers with high performance [10-12]. 
However, the self-organization process of QD array form- 
ing is difficult to manage. The random inhomogeneities 
that usually exist in heterostructure essentially affect on 
this process. They lead to inhomogeneous broadening of 
spectral lines [13] and, correspondingly, to degradation 
of laser radiation quality. 

The QD design method based on the self-organization 
effect is, in fact, an alternative to electron lithography 
whose level of development does not permit to imple- 
ment the ordered array of approximately the same QDs 
with small enough dispersion of their sizes. In the present 
paper we propose another way to design semiconductor 
laser system. This way presupposes the usage of the 
same regular microscopic areas ordered array that can be 
implemented by lithography as an active laser medium. 
This proposal is based on the following. At present tech- 
nology of GaAs mono-crystal with super high mobility 
and big length of phase coherence of charge carriers that 
achieves 10 iim and more is well-developed. Due to that 
the system of microscopic regular areas with potential 
barrier on the margin of each of them made from such 
materials can be implemented. Because of big length of 
phase coherence, carriers will take part in ballistic mo- 
tion and reflect back from area margins. The motion of 
charge carriers is similar to dynamics of billiard systems. 
Since this motion is described by Shrodinger equation 
such a billiard system can be characterized as a quantum 
billiard (QB). 

Due to finite motion of charge carrier, electron spec- 
trum of QB is discrete. At the same time it is quite dense, 
because of a QB size is much bigger than the wavelength 
of quasi-particles. This fact makes complicate the QB 
usage as an active system for the semiconductor laser, 
because dense spectrum decreases frequency stability of 
laser radiation. The frequency jumps appear easily at 
the small deviation of control parameters in the laser res- 
onator with dense frequency spectrum. As a result there 



is a problem of "rarefaction" of the spectrum: can the 
dense QB spectrum be done much sparser without chang- 
ing of QB geometrical parameters? 

A quasi-optical cavity resonator has dense and discrete 
frequency spectrum as well. The Maxwell equation de- 
scribing electromagnetic oscillations in such a resonator 
coincides with corresponding scalar Shrodinger equation 
at definite conditions. All that gives an opportunity to 
use the quasi-optical resonators as model objects to study 
spectral properties of QB. Similar modeling using mi- 
crowave resonators was done earlier for study of the phe- 
nomena that is relevant to quantum chaos [14-17]. 

In the present paper the nature of spectral lines broad- 
ening and their shifting are studied in a quasi-optical 
millimeter wave cavity resonator with random inhomo- 
geneous inside. To clarify the mechanism of broadening 
and shifting the original statistical spectral theory based 
on the separation mode technique was worked out. This 
technique was designed earlier by one of the author for 
open waveguide systems [18,19]. Using this technique we 
found out that the mechanism of spectral lines broad- 
ening was caused by intermode scattering. At that the 
largest broadening are subjected to closely set lines. And 
solitary lines keep high quality factor and intensity in the 
random inhomogeneous resonator. It results in so-called 
spectral " rarefaction" . 

The systematic experimental spectral measurements 
of numerous realizations of quasi-optical random inho- 
mogcneity cavity resonator filling were carried out for 
theory results implementation. The experimental results 
proved that the main mechanism of broadening and shift- 
ing of spectral line is inter-mode scattering caused by 
inhomogeneities resonator filling. We detected theoret- 
ically predicted effect of "rarefaction" of the resonator 
spectrum. Thus, we determined that the influence of in- 
homogeneities on the resonator spectrum has both, as 
known, negative nature (broaden lines) and a positive 
one. This positive influence is essential "rarefaction" of 
the spectrum. This fact can be important in QB usage as 
an active system in a semiconductor laser approach. Here 
we state that the high quality modes are not subjected, 
practically, to broadening and can satisfy self-excitation 
laser condition even under big inhomogeneities in an ac- 
tive resonator medium. Exactly the stable laser genera- 
tion can be implemented under one of such modes. 

Another aspect of influence of inhomogeneities on res- 
onator spectrum relates to its stochastization. In the 
paper we carried out statistical analysis of the resonator 
spectrum with a different number of inhomogeneities. We 
found out that the inter-frequency (IF) interval distribu- 
tion in the resonator spectrum is close to Poisson dis- 
tribution in a small number of inhomogeneities. This 
kind of distribution has systems with non-correlated IF 
intervals. The stochastic spectrum part appears at the 
increase of the number of inhomogeneities. The spec- 
trum is mixed, i.e. contains regular and stochastic parts. 
We estimated relationship between regular and stochastic 
spectrum parts by statistical distribution analysis of IF 
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FIG. 1: The geometry of the cyhndrical quasi-optical cavity 
resonator. S is the resonator side face, Q is its volume, H is 
the height of the cylinder, R is the radius of its base. 



intervals depending on the quantity of inhomogeneities in 
the resonator. We carried out the modeling of semicon- 
ductor laser based on QB with random inhomogeneities. 
For the purpose of this we used a quasi-optical millime- 
ter wave cavity resonator filled with random inhomo- 
geneities and with an inserted Gunn diode as an active 
element. In such a self-oscillation system we studied self- 
excitation oscillation conditions. We found out that the 
generation is unstable and multi-frequency exists in an 
empty resonator near the excitation threshold because of 
dense spectrum. Such generation can be explained by fre- 
quency jumps in dense resonator spectrum. These jumps 
disappear because of " rarefaction" of the spectrum if in- 
homogeneities are inserted into the resonator, and total 
number of generating frequencies is decreased at a defi- 
nite range of control parameters deviation. 



II. STATISTICAL SPECTRAL THEORY OF 
THE RESONATOR WITH RANDOM BULK 
INHOMOGENEITIES. 



A. Statement of the problem 

Let us consider a cylindrical quasi-optical cavity res- 
onator of radius R and height H (see FigQ. Inner vol- 
ume of the resonator, Q, is assumed to be filled with 
the material having random inhomogeneous permittiv- 
ity. We will be interested in oscillations that turn into 
transverse-electrical resonance mode (TE'-mode) pro- 
vided that the resonator is filled homogeneously. The 
vertical (z) component of the electrical field of this mode 
is equal to zero. 

According to Ref. |23|, electromagnetic field of the 
TE mode can be calculated trough magnetic Hertz vec- 
tor having only one non-zero component, namely, z- 
component n2(r). We assume that inhomogeneity of per- 
mittivity of the resonator infill is small. In this case to 
define II^ (r) in the inhomogeneous resonator we can use 



the approximate wave equation 

[A + k^e{r)] n,(r) = 



(1) 



where all components of the Hertz vector, except z- 
component, are assumed of zero value taking into ac- 
count small inhomogeneity. In Eq.J^I, A is the three- 
dimensional (3D) Laplacian, £(r) = £o + Se{r) + ia is the 
complex permittivity whose imaginary part a takes phe- 
nomenologically into account ohmic loss in the system; 
the function 5e{r) describes random space fiuctuations of 
the permittivity around its average value Eq, k = lu/c is 
wavenumber. 

In this paper we are interested in the analogy be- 
tween electromagnetic resonators and solid-state micro- 
objects, namely, QDs and QBs. In the case of classical 
resonance system, the problem of excitation by a given 
point monochromatic source is governed by the equation 
Eq. which should be supplied with S-teim in the right- 
hand side. The equation thus obtained coincides in form 
with the equation for Green function of quantum parti- 
cles moving in dissipative medium and being subjected 
to inhomogeneous scalar potential. Having this in mind 
we will carry out all further analysis relevant to the elec- 
tromagnetic resonator in terms of the dynamic equation 
for particles in the quantum dot of cylindrical shape. 

The equation for Green function of quantum particles 
moving in dissipative medium under the influence of ran- 
dom static potential V{r) has the form 



[A -f fc2 _ i/Td - V{r)] G(r, r') = ^(r - r' 



(2) 



Here is the dissipative attenuation time whose inverse 
value has the same physical meaning as the imaginary 
part of function £(r) in Eq.l^ taken with minus sign. 
The potential V{r) in the case of electromagnetic sys- 
tem is given by V{r) = —k'^Se{r). Boundary conditions 
for the solution to Eq.© result from the requirement of 
vanishing tangential components of the electrical field of 
TE mode on the resonator interface. At side boundary 
5, the Neuman condition 



5G(r,r') 



dr 



= 0, 



should be met whereas at end surfaces z 
Dirichlet condition 



G(r,r' 



z=±H/2 



= 



(3a) 
±H/2 the 

(3b) 



should be satisfied. 

For studying oscillations spectrum of the resonator 
with random inhomogeneities, the poles of Green func- 
tion averaged over realizations of the potential V^(r) from 
Eq.® should be determined. This function can be found, 
for example, from Dyson equation. Yet there by now do 
not exist effectual methods for solving this equation in 
the case of confined multi-dimensional systems. In this 
paper, we apply for this purpose the original calculation 
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technique which rehes on precise separation of modes in 
an arbitrary confined system, including the disordered 
one. The method of mode separation was previously 
developed for solving transport problems in disordered 
2D open systems and then modified fo r sy stems of 
waveguide geometry in three dimensions [I^, UM, ■ 
low we set forth adaptation of this technique for systems 
of closed geometry, in particular, for cavity resonators 
and quantum dots. 

Let us at first turn to mode representation of the equa- 
tion Eq.Q using some set of basis functions. The most 
appropriate for our purposes seems to be the whole set 
of eigenfunctions of the Laplace operator. For the cylin- 
drical resonator shown in FigQ these functions can be 
factorized to the form 



|r,M> = \r,(p;l,n) \z,q) 



(4) 



where r = (r, ip, z) is the radius-vector in cylindrical co- 
ordinates, /X = {l,n,q) is the vectorial mode index conju- 
gate to that vector. Normalized eigenfunctions of "trans- 
verse" part of the Laplacian, which obey boundary con- 
ditions Ea. lp^ajl . are given by 

|r, ^; /, n) - Cm/ (V^R) J|„| (7^'"' V/i?) e™'^ (5a) 
1 = 1,2,..., n = 0,±l,±2,... , , 
where the coefficient Cin has the form 



Cln — 



7/ 



(l»l) 



7; 



(l«l) 



-,1/2 



(5b) 



Ji 



|n| 17, 



,(l"l) 



The set of coefficients 7}'"" in Eqs.jSJ consists of positive 
zeros of the function J^^^^ (t) which are numbered by in- 
dex I in ascending order. The eigenvalues corresponding 

to functions Ea. (|^ are equal to = — (7;""'''/ ^) • 

Basis functions of the "longitudinal" Laplacian, d^/dz^, 
which meet boundary conditions Ea. lj3b|l . are given by 



9 = 1,2, 



sm 



z 
H 



nq 



(6) 



the corresponding eigenvalues being equal to (nq/H) . 
In the basis of functions Q , Eq. ^ takes the form 

(fc^ - kI,- i/Td- V^) Gfj^fj.' - ^ Ufj,^Gi,fj,' = (5^^' . 

(7) 

Here, G^^' is the Green function in mode representation, 
the parameter 



.2 _ 7, 



R 



(8) 



is the unperturbed "energy" of the mode fi (the eigen- 
value of 3D Laplace operator), functions Ufj^^, are mode 
matrix elements of the random potential. 



U^^= / dr{r;ti\V(r)\r;,y) 



(9) 



Attention should be drawn to the fact that the in- 
tramode, i. e. diagonal in mode indices, matrix element 
U/j,^ = Vfj. is separated in Eq.© from other terms of 
the sum where thus only matrix elements correspond- 
ing to intermode scattering are left. It was shown in 
Ref. that such a separation of intra- and intermode 
effective potentials provides mathematical correctness of 
the derivation of closed equations for the diagonal com- 
ponents of Green matrix ||G^pi'|| and through those com- 
ponents for that matrix integrally. 



B. Separation of the modes 

To solve the infinite set of coupled equations Eq. O is 
not less intricate problem than direct solution of multi- 
dimensional differential equation Eq.®. The problem 
would be resolved elementary if the resonator modes al- 
lowed for their strict separation. Normally, modes are 
easily separable if the resonator has no inhomogeneities. 
It will be shown below, based on the technics developed 
m Refs. dill 113, that in fact they can be separated 
even in the case of arbitrarily inhomogeneous resonator. 
But the cost of this separation in the general case is the 
appearance in equations for each of the modes of the 
effective potential, known as T-matrix in quantum the- 
ory of scattering [23 , whose functional structure is much 
more involved than that of the initial potential V{r). 

As a starting point for mode separation we introduce 
unperturbed (or trial) mode propagator by omitting 
in Eq.jTI) all intermode potentials W^^, 

G^^^ = {e-Kl-t/Td-V^)-' (10) 

The term "unperturbed" will thus be used hereupon with 
respect to intermode potentials, intramode ones being 
accounted for precisely. 

By substituting /x' = /x in Eq. |(7J) we obtain linear non- 
uniform connection of intramode propagator G^^ with 
all intermode Green functions having the particular right- 
hand mode index fi. 



(11) 



Assuming then fjb' ^ jj, and performing some necessary 
re-labellings of mode indices we can reduce Eq. {7)) to the 
form 

G[,^' Gufj, — y ^ Uvv'Gi,' = lAi,fj,G (12) 
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The latter system of interconnected equations can be 
solved with respect to intermode elements of the Green 
matrix. In the operator form the solution is given by 

G^^ = (l - i?) ' itPf^Gf,^ , (13) 

where the linear operator i? = &^^U of intermode scat- 
tering is introduced, which acts in mode subspace 
consisting of the whole set of mode indices but the index 
/i.; is the projection operator whose action reduces to 
the assignment of the value /i. to the nearest mode index 
of any adjacent operator, no matter standing to the left 
or to the right of P^. Operators and U are specified 
on by matrix elements 

{v\G^^^ \v') ^ G^^H^^, , (14a) 
{v\U\i^') ^U,,,,, . (14b) 

Correspondingly, matrix elements of the operator R are 
given by 

{v\R\u')^G'-^'>U^^> . (15) 

By substituting intermode propagators in the form 
Ea. (|13|l into the relationship Ea. ljllll we obtain the 
following rigorous expression for intramode propaga- 
tor G^^, 

G^^ = {k' -nl~ilTi-V^-r^y^ . (16) 

Here 

r^ = P^ri'(l-i?)"'i?P^ (17) 

is the portion of the mode eigen-energy which relates 
to the intermode scattering. 

It should be noted that for determining the disordered 
resonator spectrum the poles of solely diagonal elements 
of the Green matrix suffice to be found. Just these ele- 
ments determine, in view of relationship Ea. (|13|l . all ma- 
jor analytical properties of the whole Green function from 
Eq.©. Below we will analyze the cavity-resonator spec- 
trum in frames of relatively simple statistical model of 
random potential V(r). 

C. Statistical analysis of the resonator spectrum 

We assume the potential V(v) to have zero mean value, 
(y(r)) = , and the binary correlation function 

(y(r)y(r')) = DW{v - r') . (18) 

Bearing in mind the forthcoming numerical analysis we 
will take the function W[v) in the form of Gaussian ex- 
ponent, viz. W[y) = cxp(— r^/2r^) , where stands for 



the correlation radius. In the case of electromagnetic res- 
onator the normalization constant D in Ea. H18|l is given 
hy D = k'^a^, where = (fc^(r)) is the variance of 
permittivity fluctuations. 

The pair of selected statistical parameters, namely, the 
average random potential and its binary correlation func- 
tion, are sufficient for making detailed analysis of the 
system in study if the function fc(r) is the Gaussian- 
distributed random variable. Yet these two parameters 
suffice for carrying out the asymptotically correct anal- 
ysis even in the case where statistics of the fluctuations 
is markedly non-Gaussian, provided that the potential 
V{r) is in a certain sense small. As it is conventional in 
condensed matter physics, we will regard the potential 
to be small and the resulting scattering, correspondingly, 
weak provided that the scattering rate calculated in Born 
approximation is small as compared to the unperturbed 
quasiparticle energy, k^ in our case. The smallness of one- 
fold scattering probabihty enables one to regard the po- 
tential V{r) with parametric accuracy as Gaussian ran- 
dom process, whatever in fact its distribution is p^ . 

Consider now the self-energy operator of the mode fj, 
which consists, in accord with Ec ljltil) . of two terms, 

^^ = V^ + T^ . (19) 

The first term in the right-hand side of this formula van- 
ishes when being averaged whereas the second term does 
not. Its average value can be easily calculated under the 
presumption of weak scattering. The strength of the in- 
termode scattering is estimated by the norm of the oper- 
ator R entering Ea. p7|l . Assuming this norm to be small 
as compared to unity and keeping only two main terms in 
the expansion of the inverse operator in Eq. H17|l we come 
to the necessity of performing the averaging not of the 
exact operator potential 7^ but rather of its relatively 
compact limiting value 

« P^UG^''^UP^ = U^uG^;:'^U^^ . (20) 

When calculating the quantity (^), one can neglect 
with parametric accuracy the correlation between in- 
tramode and intermode potentials. This enables us to av- 
erage intramode propagator G'lf ^ and its envelopes con- 
sisting of intermode potentials U^^u independently. To 

average the function G'i^'' it is worthwhile to present it, 
at first, in the integral form, 

oo 

Gi^) = j dt exp [-i {k^ -kI^ i/Td - V^) t] . (21) 



Then the averaging of the integrand in Ea. (|21|l with 
the aid of continual integration with Gaussian functional 
weight yields 
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~i (1 - Kl/k^ - i/k^Td) t - -n'^L^ir,) 



exp 



2k'^a^L^{r^) 

I 



1 - $ 



(22) 



Here $(^) is the probability integral 0|, L^{rc) is the the case of Gaussian correlation function M^(r) we readily 
dimensionless correlator of intramode potential Vv In obtain 

I 



i.(rc) = (V^V,) = / / dvdT' (r, I. |r, v) exp [-(r - r' f /2rl] (r', |r', 



1 1 



—Cf^^^ j j dsds' exp 







tt2 



siv? [nqn^ s) siv? {nqn^ s') 



1 1 



X j jtt'dtdt'jf^J^lyk)jf,^J^lyk')fd^e^p 





i?2 

ttU'^ + t'^ - 2tt' COS V?) 

2ri ' 



(23) 



r 



Although the expression Ea. H22l) is formally exact, it is 
not quite convenient for the analysis in view of its bulky 
structure. Asymptotical calculations at large and small 
values of the probability integral argument permit us to 
make use of much simpler interpolation expression, 



(24) 



which is close in form to the initial unperturbed Green 
function Ea. (|10|) . Here we have used the notation 
l/r* = 1 /r^ + fc2 -y/ [2l'K)a'^Ly(r^ for the effective scat- 
tering frequency. The latter includes both the initial dis- 
sipative term l/r^ and the addendum originating from 
wave scattering by random inhomogeneities. 

Interpolation of the average trial Green function by 



J 



the expression Ea. (|24|l permits us to interpret the quan- 
tity k^ -y/ (2 1 -K^a"^ L ^(r c) as the dephasing frequency of the 
mode state which is related to scattering by the inho- 
mogeneities in the resonator. 

To make further comparison of theoretical results with 
experimental data we will consider below the specific lim- 
iting case where the inequalities are fulfilled 



rc<H<t:R 



(25) 



These restrictions enable us to calculate the integrals in 
Eq. (|23() asymptotically, thus resulting in the following es- 
timate for the parameter L^irc)'- Li,{rc) {rd R) ■ The 
correlator of intermode potentials in Eq. H20() can be rep- 
resented as 



1 1 



dsds' exp 







l\2 



2^ 



sin(7rg„^s) sin(7rq„^s') sin(7r(7„^s) sin(7rg„^s') 



1 1 



X / / tt'dtdt'J^„J^iyk)J,,J^^^^^^^^^ 







(j) <j> d(pd(p' exp |— i(n^ — n^){tf — tf') — ^-^ [t^ + i'^ - 2tt' cos((p — (^')] | 



(26) 
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Asymptotic calculation of the integrals over ip and ip' where the factor A^^(rc) is given by 
results in the following formula 



(27) 



1 1 





1 1 

X y" y VWdtdt' exp 





If 

2r 



tt2 



sin (nqn^s) sin [nqn^s') sin (nqn^s) sin (7rq„^s') 



2 1 \2 



X J, 



v,(7^'^0^«.,(7^'^t')J,„,(.^|)t)J,„,(,^|)t' 



(28) 



r 



With the results (|24|l and H27|l . separation of real and 
imaginary parts of the average T-matrix, (7^) — Afc^ + 

i/Tjf'^\ results in the following expressions for shifting 
and broadening the fi-th resonant level, 



ich) 



(29a) 



(29b) 



by its interactions with the adjacent resonances, whose 
position on the frequency axis are confined to the region 
limited by order equality ^ I/t*. In the case 

where only one resonance level with, say, frequency Wi, 
proves to fall into the above indicated interval around the 
/x-th resonance (we assume the condition ard R ^ 1 of 
weak intramode scattering to be met) its contribution to 
the shift and the width of the /x-th resonance is estimated 
as 



(31a) 



The factor A^i^(rc), as can be seen from Ea. (|28|l . is a 
real-valued quantity, its absolute value being estimated 
as rc/R <C 1. Although the sign of this factor cannot 
be uniquely identified in the general case since it con- 
tains the dependence on specific indices of modes between 
which the scattering is carried out, numerical analysis 
shows that A^j^(rc) > 0, which corresponds to the broad- 
ening of spectral lines. One can adequately estimate both 
the shift and the broadening of /x-th resonant level by 
substituting the function (gL^'') in the interpolated form 
Ea. (|24|l . instead of its exact expression Ea. H22|) . into the 
right-hand sides of formulas Ea. l29() . The result in this 
case reduces to 



tch) 



I 4 2 

— k a — 



(30a) 



(30b) 

The structure of the summands in Eas. (|30|) indicates 
unambiguously that both the shift and the broadening 
of each given resonance {k"^ ~ fc^) are mainly provided 



(31b) 



Here 



and 



are the cyclic spectral frequency of 



empty resonator and that of the resonator filled with 
random inhomogeneities, respectively, ^oj^q a-nd 5uj^ are 
relative widths of their spectral lines, the parameter 
A ^ Tc/R- If several resonances fall into the indicated 
vicinity of the /i-th resonance simultaneously, they con- 
tribute additively to both the level position and width. 

The above presented theory shows that the filling of a 
quasioptical cavity resonator with randomly distributed 
inhomogeneities results in both the random shift of its 
resonance lines and the increase of their linewidths. It 
should be noted that formulas Eas. H29|) and H30|) were ob- 
tained with the use of some procedure of averaging over 
the ensemble of macroscopically identical resonance sys- 
tems whose microscopic realizations are different. What 
this formally means is that the statistical theory can give 
only qualitative predictions for the experiment. The sit- 
uation seems to be quite similar to that one faces when 
studying problems of wave scattering by randomly rough 
surfaces jjj . Yet the problem of oscillations in the quasi- 
optical cavity resonator filled with bulk random inhomo- 
geneities differs substantially from the problem of wave 
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scattering by randomly rough surface. The point is that 
oscillations in cavity resonators become established as a 
result of multiple transmissions of waves (with multiplic- 
ity of the order of Q) through the random infill of the 
resonator. With this consideration in mind one might 
expect that under the restriction Tc <C i? the characteris- 
tics of resonance lines should be effectively self-averaged 
p^ l. If this is the case, the agreement between theory 
and experiment may be achieved either in the course of 
small number of trials or even for one particular real- 
ization of the random system. The confirmation or the 
denial of this assumption may be achieved through the 
comparison of theoretical predictions and the data of cor- 
responding measurements. 



D. Numerical simulation 



To compare the theory and experiment we calculate 
the spectrum of TE modes in the quasi-optical cylinder 
cavity millimeter wave resonator filled with random di- 
electric inhomogeneities using our theoretical results. At 
the resonator spectrum calculation the resonator param- 
eters were taken according to the experiment conditions. 
The spectrum was calculated by solving the excitation 
task with pointed external dipole. Fig. 2 shows the influ- 
ence of random bulk inhomogeneities on the spectrum at 
the different values of the parameter a calculated using 
Eq.(30). As an example we selected the spectrum in the 
frequency interval of 36-37 GHz. 

It is shown that the spectrum essentially changed at 
the presence of dielectric inhomogeneities in the res- 
onator. The magnitude of shifting and broadening of 
spectral lines are caused by resonator filling by inhomo- 
geneities. The solitary spectral lines even at presence 
of inhomogeneities keep high quality factor, and their 
broadening is enough small (for example the spectral 
lines in the range of 36.2-36.6 GHz). At the same time 
the shifting and broadening of adjacent lines are signifi- 
cant (for example, the spectral lines near of 36.8 GHz). 
This results that the number of excited lines with high 
quality factor are essentially decreased if the number of 
inhomogeneities are getting bigger. The reduction of the 
number of high quality lines we can explain as "rarefac- 
tion" of the resonator spectrum caused by unequal scat- 
tering conditions of different modes on inhomogeneities. 
Fig. 3 shows the amplitude-frequency dependence of the 
resonator spectrum built as the frequency dependence of 
Green function module at different magnitude of param- 
eter (T. So, at the increase of the number of inhomo- 
geneities in the resonator the adjacent resonances over- 
lapped and the quality factor of combined resonances is 
getting smaller (for example, spectral lines near by 36.8 
GHz). 
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FIG. 2: The dependence of relative shifting P — Ak^^/k^ = 
(a,c) of spectral lines of the resonator on frequency / 
caused by bulk inhomogeneities and their relative broadening 
W = T^/rjf'''' — Qd/Qii (b,d) at the permittivity dispersion 
values of inhomogeneities: a = 0.02 (a,b) and a = 0.05 (c,d); 
Vc = 0.3 cm. 



III. EXPERIMENT AND DISCUSSION 

The main goal of our experiment is the verification of 
the numerical simulation results of shifting and broad- 
ening of spectral lines caused by inhomogeneities and 
the possibility of resonator spectrum " rarefaction" . The 
another goal is as follows. In the paper [3] we de- 
tected strong spectrum stochastization caused by inho- 
mogeneities such as anisotropic sapphire particles with 
the dimensions of order of operating wavelength placed 
into the cavity resonator. The spectrum has mixed state 
due to both regular and chaotic spectrum part exist. In 
contrast to [3], in the present paper we study the influ- 
ence on the resonator spectrum relatively small isotropic 
inhomogeneities. Such inhomogeneities can be made by 
styrofoam particles with the average value of the permit- 
tivity closed to one and with small dielectric loss angle. 
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FIG. 3: The amplitude-frequency dependence of the empty 
resonator spectrum (a) and the spectrum of resonator with 
inhomogeneities at a = 0.02 (b) and a = 0.05 (c); Vc = 0.3 
cm. The amphtude is normahzed by maximal amplitude of 
resonances for the empty resonator in the considered wave 
range. 
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FIG. 4: The quasi-optical cylinder cavity millimeter wave 
resonator filled with inhomogeneities, 1,6 are input/output 
waveguides, 2,5 are holes coupling with waveguides, 3 is the 
resonator body, 4 is styrofoam particles, D is the resonator 
diameter, H is its hight; D = 130 mm, H = 14 mm. 



Thus, there is a question: is it possible the spectrum 
chaotization here? 



A. Experiment technique 

The study of the influence of inhomogeneities on the 
cavity resonator spectrum we carried out at frequency 
range of 32-37 GHz. We took a quasi-optical cylinder 
millimeter wave resonator random filled with styrofoam 
particles (Fig. 4). The styrofoam particles have the real 
part of permittivity is about one, close to the open-air 
permittivity and small dielectric loss, e — 1.04 -I- ilO^'*. 

To study the influence of small dielectric particles on 
the resonator spectrum it is necessary to provide high 
quality factor for the resonator oscillations without inho- 
mogeneities. For that we excited TE mode. The mag- 



netic field vector of this mode is directed along the res- 
onator z-axis, and microwave currents do not cross the 
interface between the flat resonator face and cylinder sur- 
face. Owing to that the empty resonator has high qual- 
ity factor up to 2 X 10^ . To excite the selected mode 
we used a waveguide diffraction antenna. It is the circu- 
lar hole with the diameter of 2 mm in a thin diaphragm 
with the thickness of 0.1 mm closing the input waveg- 
uide. The diaphragm surface is flush-mounted with the 
side-cut cylinder surface of the resonator. The same an- 
tenna is used to receive the oscillations on the opposite 
side of the cylinder surface. 

The resonator spectrum was detected using "on pass" 
regime in 32 - 37 GHz wave range by the wide frequency 
standing wave ratio meter. Measurement process was 
automated. Signal from the measurement device using 
analog-digital conversion goes into computer. The fur- 
ther signal processing (the determination of the spectral 
line intensity, its quality factor, and frequency) was done 
by special software GUI application. It gives resource 
in a short space of time to handle measurement data 
for a huge number of realizations of the random inhomo- 
geneities distribution in the resonator. Owing to that the 
frequency and quality factor measurement accuracy for 
spectral lines do not exceed 0.1% and 1%, respectively. 
The styrofoam particles that are used as inhomogeneities 
have the size about of 2-3 mm. The space distribution 
of them was arbitrary for each realization. The spectral 
characteristics were measured depending on the number 
of these inhomogeneities. 



B. Shifting and broadening of spectral lines. Effect 
of spectrum "rarefaction" 

The spectrum of the empty resonator is dense enough 
and consists of 84 narrow spectral lines in the range of 
32-37 GHz (Fig. 5). Each line was identified according 
to mode indexes of corresponding to its eigen resonator 
oscillation. We found out the oscillations with high qual- 
ity factor (Q) of order of 10"* and higher have small az- 
imuth indexes (n 1) and high radial indexes (/ >> 1) 
. For the base oscillation mode (n >> 1,^ = 1) that 
has the field distribution concentrated inside of the res- 
onator side-cut (whispering-gallery oscillations [5]), Q is 
about 2 • lO'^. The spectrum of the resonator with in- 
homogeneities is essentially different from the spectrum 
of the empty one. By the increase of the number of in- 
homogeneities spectral lines are getting wider, and Q is 
getting smaller. The intensity of broadening lines is de- 
creased as well. The spectrum keeps only few lines with 
enough high intercity and Q value. For such lines the 
Q value is close to the Q-factor of the spectral lines in 
the empty resonator. So, the number of spectral lines 
with enough high intensity and with Q > 10^ is 52 in the 
empty resonator (Fig. 5a), at filling the resonator by sty- 
rofoam particles the number of such lines is 10 (Fig. 5b), 
and at filling the resonator by pressed styrofoam particles 
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FIG. 5: The spectrum of the empty resonator (84 spectral 
lines) (a); the resonator filled with styrofoam particles (77 
spectral linos) (b); the resonator filled with pressed styro- 
foam particles (57 spectral lines) (c); The normalization of 
amplitude was made on the maximal amplitude value for the 
empty resonator. 



the number of lines is 3 (Fig. 5c). 

The presence in the resonator with inhomogeneities to- 
gether with broadening lines few high Q spectral lines is 
equivalent to the spectrum "rarefaction" (Fig. 5). It is 
necessary to note that the presence of lines with high Q 
indicates that the broadening is specified, mainly, by the 
inhomogeneous of the resonator filling permittivity and 
is not specified by additional dissipation loss caused by 
small dielectric loss in the styrofoam. Except broadening 
of spectral lines they have frequency shifting. This fre- 
quency shifting has both regular component and random 
one as well (Fig. 6). 

The regular shifting happens towards the low fre- 
quency direction because of the increase of average di- 
electric permittivity of inhomogeneous medium in the 
resonator. For example, it is 150 MGz for the resonator 
fully filled with styrofoam particles. 

It is significant that if the distance between resonances 
is enough small they can both come close and depart 
from each other because of influence of inhomogeneities. 
Their Q is essentially changed (Fig. 7, 8): it can increase 
for one of them and decrease for another. 

Both the spectral lines broadening and their shifting 
have quality agreement with described above theory and 
can be interpreted in the terms of intermode scattering. 
The maximal broadening is for the adjacent lines that is 
agreed with theory. 

We can give the following physical explanation of spec- 
tral lines broadening caused by random inhomogeneities 
based on quantum mechanics interpretation of resonator 
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FIG. 6; The dependence of the frequency shifting 5f = 
U empty - fin hom) /f empty On frequency / for the resonator 
filled by inhomogeneities. The dotted line is the value of shift- 
ing regular component. 
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FIG. 7: The relative Q-factor deviation 5Q = {Qempty — 
QinhoTn)/Qe,npty dcpcnds upon frequency / for the resonator 
filled by inhomogeneities. For major part of resonances the 
relative Q-factor deviation is more than zero, i.e., their Q is 
less than for the empty resonator. The relative Q-factor devi- 
ation is less than zero for several nearest adjacent resonances 
that have overlapped spectral lines (in the case of the empty 
resonator) and split into separate resonances at the presence 
of inhomogeneities, i.e., the spectral line repulsion takes place 
here (see Fig.8). 



oscillations. In empty resonator the lifetime of a reso- 
nant mode in a given state (describing by an eigen reso- 
nant frequency) and its spectral line width define by dis- 
sipation loss in the resonator body. The uncertainty of 
mode energy within the spectral line width is caused by 
this loss. The intermode interaction appears at insert- 
ing inhomogeneities into the resonator. Owing to that 
the uncertainty of mode energy and, correspondingly, the 
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FIG. 8: The adjacent resonances with indexes 
n = 34, 1 = 3; n = 4, 1 = 15; in the empty resonator (a) 
and in the resonator filled by inhomogeneities (b). The 
normalization of the spectral lines amplitude was made on 
the maximal spectral line amplitude for the empty resonator. 



spectral line width increase due to the mode transition 
between nearest states. The mode energy loss due to dis- 
sipation loss in the resonator body and transition to the 
neighboring state are additive to the spectral line width. 



C. Statistical analysis of IF intervals 

If the number of inhomogeneities increases the spec- 
trum of the quasi-optical resonator possesses stochastic 
character that is visualized in distribution of IF inter- 
vals. In order to define the relation between regular and 
random spectral components the comparison of IF inter- 
vals distribution obtained experimentally with different 
theoretical distributions based on a priori data about sta- 
tistical process is usually used. In particular, we use the 
Brody function determined the distribution of IF inter- 
vals probability -Pb(s) that is given by 



(32) 



where s — (tOn — w„_i)p(a;„), tOn is the spectral line 
frequency, is the spectral density - the superposi- 

tion of regular and random motion density, /3 is the mea- 
sure of stochastic motion, constants A and B are defined 
from the condition of standardization: A = {1 + (3)3, , 
B = T^+P{2 -h /3)(1 + r{z) is the Gamma function. 
At /? — > IF intervals in the spectrum are not correlated 
and can be described by the Poisson distribution, and at 
/? ^ 1 we have Wigner distribution, when the repulsion 
eff'ect of spectral lines exists, that is the probability of 
closest to zero inter-frequency interval is equal zero as 
well. 



If stochastization measure in the spectrum is relatively 
small and the distribution of IF intervals is close to 
the Poisson one we can use Berry-Robnik distribution 
Pbr{s) [26] that is given by 



P, 



(s) = p^e-P'erfc (1 - p) + 
^{l-pfs + 2p\ (1 - p) e-P^-f (I-'')' 



(33) 



where erfc(a;) = J exp(— t^)dt. 

X 

p = 1 is relative phase volume occupied with regular 
trajectory in mixed systems. The limit p ^ 1 corre- 
sponds to regular system; p — > is completely chaotic 
system. 1 — p is relative phase volume occupied with 
chaotic motion. 

The experimental data show that in the empty res- 
onator the IF intervals distribution is close enough to the 
Berry-Robnik distribution with p = 1 (Poisson distribu- 
tion, P(s) exp(— s))(Fig.9a). By increase of the num- 
ber of inhomogeneities (styrofoam particles) the function 
P(s) is reduced and it has maximum at small s. The 
presence of the maximum P{s) indicates the appearance 
of the repulsion effect of spectral lines, and the random 
component is increased. We found out that at full filling 
the resonator by styrofoam particles p = 0.4 (Fig. 9b). 

As we can see in Fig. 9 the function P{s) has maximum 
at small s that is not described by the Berry-Robnik de- 
pendence [26] . The presence of such a maximum can be 
explained by Chaos- Assisted Tunneling (CAT) [27] in the 
resonator with random inhomogeneities. The proposed in 
[27] IF distribution function gives the possibility to de- 
scribe the observed distribution at the parameter p = 0.4 
corresponding to the classical dynaniicrs in the system and 
the parameter = 0.1 describing the tunneling between 
different states of modes. 

We calculated the spectral rigidity for the resonator 
with inhomogeneities. The spectral rigidity A3(L) is an 
integral characteristic of degree of spectral lines ordering 
for frequency distances that is much more than inter- 
frequency interval. It is given by [15] 



x+L 



A3{x,L) = ^mm j [n{e) - As - B]'^de 



(34) 



where L is the interval length on which the function 
A3(a;, L) is determined. The function n(e) is built as 
follows [15]. For the sequence of the frequencies lo^ nor- 
malized on unit density (a;„ = co^_i + Sn), we introduce a 
staircase function n{e) equals the number of frequencies 
with LO^ < s. 

The function n{£) has a staircase view with aver- 
age unit tilt. The function A3{x,L) is determined as 
the minimum of quadratic deviation n(e) in the inter- 
val (x, x + L) from the straight line. The meaning of 
the averaged in x spectral rigidity (Eq.(34)) (A3(a;, L))^ 
depends only on L and is denoted as A3(L). 
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FIG. 9: The distribution of IF intervals, P(s). Figure 9a 

is for the empty resonator spectrum. Sohd hue is for the 
Berry- Robriik distribution at p ^ 1 or the Brody distribution 
at /3 ^ 0. Figure 9b is for the spectrum of the resonator 
fully filled with styrofoam particles. The curve 1 is for Brody 
distribution at /3 = 0.1; the curve 2 is for Berry- Robnik dis- 
tribution at p = 0.4; the curve 3 is for Podolskiy-Narimanov 
distribution at p = 0.4 and i/ = 0.1 [27]. 



The obtained curve A3{L for the resonator with ran- 
dom inhomogeneities is shown in Fig. 10. This curve (2) 
is placed between the spectral rigidity for the Poisson dis- 
tribution (curve 1, (A3(L) = L/15))) and curve 3 that is 
corresponded to the spectral rigidity for Gauss orthogo- 
nal ensemble (GOA) A^iL) = l/n^lnL - 0.00687 [15]. 
The latter is observed at the modeling of quantum chaos 
in microwave cavity unstable resonators similar to Sinai 
and Bunimovich billiards [14]. 




O.OJ ' ' ' 1 

5 10 15 20 
L 

FIG. 10: The spectral rigidity for the resonator filled with 
styrofoam particles (curve 2). The straight line is for the 
spectral rigidity of Poisson distribution (curve 1), curve 3 is 
for GOA. 



IV. ACTIVE RESONATOR WITH RANDOM 
INHOMOGENEITIES 

Wc studied also an active quasi-optical resonator with 
random inhomogeneities. In comparison with the con- 
sidered above passive resonator the active one is a self- 
sustained oscillation system where the presence of ran- 
dom inhomogeneities affects on the excitation of res- 
onator oscillations. For that wc used the same resonator 
as mentioned above (Fig. 4) with pointed Gunn diode 
inside. The microwave electrical field for TE mode was 
directed along diode's axis. The diode DC power supply 
was implemented through a filter as a quarter-wave mi- 
crowave isolator; owing to spurious microwave radiation 
was prevented. 

Quasi-optical resonator with Gunn diode is an active 
oscillator with distributed parameters. Near threshold of 
excitation in such an oscillator with the empty resonator 
was detected unstable multi-frequency generation. Such 
kind of generation we can explain by frequency jumps 
between adjacent spectral lines with high quality factor. 
If excitation threshold was highly exceeded, as a result 
of frequency competition, corresponding to these spec- 
tral lines, mono-frequency generation occurs (Fig. lib). 
The active oscillator selects "itself the only frequency 
to provide maximal regeneration factor [28]. 

The random inhomogeneities lead to effective rarefac- 
tion of the resonator spectrum. The number of adja- 
cent spectral lines with high quality factor is reduced, 
and, as a result, multi-frequency generation disappears. 
At small exceeding of the generation threshold the noise 
generation appears (Fig. 11c) in the resonator with inho- 
mogeneities. At big exceeding of the generation threshold 
the stable mono- frequency is observed (Fig. lid). Owing 
to inhomogeneities mono-frequency generation possesses 
greater frequency stability, and the number of generat- 
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FIG. 11: Oscillogram of the 36 GHz generation obtained by 
a millimeter wave spectrum analyzer, (a) and (b) are for the 
empty resonator; multi-frequency generation near the gener- 
ation threshold (a) and mono-frequency generation much far 

from the generation threshold (b). (c) and (d) arc for the 
resonator filled with inhomogcncitics; chaotic generation near 
the generation threshold (c) and mono-frequency stable gen- 
eration much fax from the generation threshold (d). 




FIG. 12: The generating frequencies / depending on the 
number of inhomogcncitics M. The empty resonator (M = 
0), the resonator quarter-filled with inhomogeneities (M = 1), 
the resonator half-filled with inhomogeneities (M = 2), the 
resonator fully filled with inhomogeneities (M = 3). 



ing frequencies is reduced in the range of diode negative 
resistance. Fig. 12 shows generating frequencies for the 
empty resonator and with random bulk inhomogeneities 
one. 



V. CONCLUSION 

The statistical spectral theory of quasi-optical cavity 
resonator filled with random dielectric inhomogeneities 



was developed in the present paper. We showed that the 
presence of inhomogcncitics leads to the broadening and 
shifting of spectral lines. It is found out that the nature 
of broadening and shifting of spectral lines is relevant 
to the intermode scattering. The scattering effect for 
the given spectral line essentially depends on frequency 
distance between of it and adjacent ones and is sharply 
decreased for bigger distances. Under the influence of 
ransom inhomogeneities original spectrum modification 
occurs that can be interpreted as spectrum rarefaction. 
The spectrum is rarefied because of solitary spectral lines 
are not practically subjected to the infiuence of inhomo- 
geneities. The quality factor of such lines and, corre- 
spondingly, their intensity stay high at the resonator ex- 
citation. The intensity of adjacent lines broadened under 
the influence of inhomogeneities is essentially reduced. 
Owing to that at the great number of inhomogeneities 
the resonator spectrum is "rarefied" - few solitary high 
quality factor spectral lines prevail in the spectrum. The- 
oretical prediction of broadening and shifting of resonator 
spectral lines and spectrum "rarefaction" are subject to 
experimental check. 

For that purpose we studied experimentally in 8- 
millimeter wave range the spectrum of the quasi-optical 
cavity resonator filled with random small-scattered bulk 

inhomogeneities. These inhomogcncitics were styrofoam 
particles with smaller size than operating wavelength. It 
is found out that such inhomogeneities lead to broaden- 
ing and shifting the spectral lines. As experiment showed 
the maximum of their influence was on frequency adja- 
cent spectral lines. The solitary lines, according to our 
theory, were subjected to this influence in much smaller 
degree. We detected the effect of stochastic spectrum 
"rarefaction". They prove that the main mechanism of 
broadening and shifting of spectral lines is relevant to 
inhomogeneities intermode scattering. We studied also 
chaotic properties of oscillations in our resonator. It is 
found out that the empty resonator has IF intervals dis- 
tribution similar to Poisson distribution that is typical to 
the spectrum with non-correlated IF intervals. Even at 
small number of inhomogeneities the resonator spectrum 
has random part that increases in proportion to the num- 
ber of inhomogeneities in the resonator. At that IF in- 
tervals distribution in random inhomogeneous resonator 
is described by Brody and Berry-Robnik distributions of 
IF intervals. 

We obtained the results concerning the influence of 
bulk inhomogeneities on the process of generation in a 
self-oscillatory system. The self-oscillatory system was a 
quasi-optical millimeter wave cavity resonator containing 
inhomogeneities with an active element as Gunn diode. 
We detected that inhomogeneities leads to the essential 
"rarefaction" of the spectrum and creates conditions for 
monochromatic stable generation in self-oscillatory sys- 
tem. The inhomogeneous quasi-optical cavity millimeter 
wave resonator (passive and active) can serve as a model 
of semiconductor quantum billiard. Based on our results 
we suggest using such billiards with spectrum rarefled by 
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random inhomogeneities as an active system of semicon- ductor laser. 
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